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Abstract 

We study two-stream instabilities in a non-equilibrium system in which a stream of electrons is 
injected into doped graphene. As with equivalent non-equilibrium parabolic band systems, we find 
that the graphene systems can support unstable charge-density waves whose amplitudes grow with 
time. We determine the range of wavevector q that are unstable, and their growth rates. We find 
no instability for waves with wavevectors parallel or perpendicular to the direction of the injected 
carriers. We find that, within the small wavevector approximation, the angle between q and the 
direction of the injected electrons that maximizes the growth rate increases with increasing |q|. We 
compare the range and strength of the instability in graphene to that of two and three dimensional 
parabolic band systems. 
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I. INTRODUCTION 


Monolayer graphene consists of a single monolayer of carbon atoms arranged in a honey¬ 
comb lattice. There has been a tremendous amount of interest in graphene, both theoretical 
and experimental, since the initial fabrication by Novoselov et al .was reported- Most inves¬ 
tigations on the electronic properties of graphene have concentrated on the linear response 
regime, in which the system is slightly perturbed from its equilibrium state (see, e.g., Refs. 

and J). In this paper, we theoretically investigate a non-equilibrium situation, in which 
a stream of carriers is injected into doped graphene sample. We show that this results in 
an instability in the collective modes of the system, which is analogous to the two-stream 
instability in classical plasmas. 

Under certain conditions in classical plasmas, when there are counter-streaming charged 
particles, some of the collective modes (he., charge density waves or plasmons) of the plasma 
become unstable, in the sense that they initially grow exponentially. This phenomenon, 
called the two-stream instability,- can qualitatively be understood by considering Landau 
damping, the process in which plasma waves decay in equilibrium. In equilibrium situations, 
the interaction between plasma waves and the charged particles in the plasma result in a net 
transfer of energy from the plasma waves to the individual charged particles, which leads to 
the decay of the plasma waves. However, in certain non-equilibrium situations the inverse 
can occur; the plasma waves absorb a net amount of energy from the charged particles and 
the waves grow in amplitude. 

The possibility of two-stream instabilities in solid-state systems has been studied theo¬ 
retically by several investigators. Pines and Schrieffer- considered the possibility of these 
instabilities in 3-dimensional solid-state systems in which both electron and holes are present, 
such as semimetals or small band-gap semiconductors. A static electric field would cause the 
electrons and holes to counter-stream in opposite directions and can in principle produce the 
instability. However, in practice, the strong electron-hole scattering in 3-dimensional systems 
suppresses the counter-streaming motion of the electrons and holes. Several groups have the¬ 
oretically investigated two-stream instabilities in coupled two-dimensional structures such 
as closely-spaced electron and hole doped quantum-wells to separate the oppositely charged 
carriers.- - — However, even in these systems with reduced electron-hole scattering due to 
the spatial separation between the two carrier species, the strength of electric fields neces- 
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sary to obtain a sufficient large relative drift velocities of the electron and hole populations 
would cause heating of the carriers that suppresses the instability. Another way of obtaining 
counter-streaming carriers with sufficiently large relative drift velocities is by injection of car¬ 
riers at high velocities into a doped system. However, in order to achieve high velocities, for 
parabolic-band systems, carriers must be injected at high energies. These high-energy car¬ 
riers usually scatter inclastically very quickly, typically with optical phonons, which makes 
it difficult to set a steady-state system with counter-streaming carriers. 

In systems with bands that have linear dispersions such as graphene, the velocity of 
a carrier is independent of the energy. Thus, it is possible to obtain large relative drift 
velocities without having to inject carriers with large energies or apply very large electric 
fields. Therefore, the reasons given above which tend to suppress the instability in parabolic- 
band systems do not apply in linear-band systems. In this paper, we show that it is in fact 
possible to obtain two-stream instabilities in linear-band systems such as graphene. 

The outline of this paper is as follows. In Section II, we discuss the formalism that we use 
to calculate the dispersion and growth rates of unstable collective modes (i.e., plasmons) 
both in a doped parabolic band semiconductor (for comparison) and in doped graphene, 
both with an injected stream of carriers. In Section III, we review the two-stream instability 
in two- and three-dimensional parabolic-band systems, and in Section IV, we investigate 
the two-stream instability in graphene and compare the results to parabolic band systems. 
Section V contains our discussion and conclusion. 

II. FORMALISM 

In order to obtain the dispersion relation and growth rates of the plasmons, we analyze the 
relative dielectric function e(q,uj) of extrinsic graphene in the presence of an injected stream 
of carriers. The dielectric function is defined to be the ratio of an externally induced potential 
V ext (q,co) to the total potential V tot (q,co) (the sum of V ext (q,co) and the internal potential 
due to the charge density perturbation in the system),— e(q,uj) = V ext (q,co)/V to t(q,uj). The 
plasmon modes of wave-vector q are obtained by solving for u q in the equation e(q,cu q ) = 0, 
since this condition implies the existence of a V to t (< 7 ,cu) without an external perturbation, 
indicating a self-sustaining wave. Since the time and spatial dependence of the density 
perturbation of the collective mode is n q exp (i [q ■ r — c o q t]), if Im(o;g) > 0 then the density 
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perturbation grows exponentially (at least initially, until non-linear effects become apparent). 

The dielectric function e(q, oj) can be obtained at various levels of approximation. We 
use the long-wavelength limit of the random-phase approximation,— which can be derived 
by solving the quantum-mechanical equations of motion of the charge carriers under the 
influence of a mean field electric potential, and Poisson’s equation. The long-wavelength 
limit gives results that can also be derived using classical plasma physics, together with the 
relationship between the velocity and momentum; i.e., in the case of parabolic band systems, 
v(p) = p/m , where m is the effective mass of the carriers, and, in the case of graphene, 
v(p) = v 0 p/\p\, where v 0 is the speed of the electron in graphene. 

We assume that our system is doped at a density ri\ of carriers, and there is no magnetic 
field. We also assume that the injected beam of carriers of density n 2 is peaked around a 
particular momentum p$. In the random phase approximation, the polarizability of the two 
charge carriers confined to the same space is equal to the sum of their individual polariz¬ 
abilities, and hence the dielectric function has the form 1 ® 


e(q,uj) = 1 - W(g)[IIi(q,u;) + II 2 (q,u)} , 


( 1 ) 


where IIi(qr, a;) and H 2 (q,oj) are the polarizability of the equilibrium carriers in the doped 
system and the injected carriers respectively, and V c (q) is the Fourier transform of the bare 


Coulomb interaction (statically screened by a background dielectric) with appropriate to 
the dimension of the system considered. We shall concentrate on long wavelength (small q) 
collective modes. In this limit the polarizabilities are given by 



( 2 ) 


which is the small q expansion of the random phase approximation.— Here, f(p) is the 
distribution function, with the normalization condition n — f dp f(p), and v(p) = V p E, 
where E(p) is the band kinetic energy. 

We assume that the doping carriers (subsystem 1) has an equilibrium distribution at zero 
temperature (be., j\(p) = constant for |p| < pp and f(p) = 0 for |p| > pp ) which is also 
appropriate when the temperature is much less than the Fermi temperature. Then in the 
small-g, non-zero frequency limit 



( 3 ) 


4 




where C is a constant which is dependent on the system dimension and energy dispersion 


E(p). 


For the injected carriers, we assume that the distribution is constant within a certain 
“distance” (in momentum-space) of p 0 and zero otherwise, which in the small q limit can 
be approximated by / 2 (p) = n 2 5(p — p 0 ). The form of n 2 (q, u>) is dependent on the system 
dimension and the E(p). We also assume that / 2 (p) retains its 5-function distribution over 
the length of the system. This assumption is of course an approximation, because this is a 
non-equilibrium distribution that is subject to collisions that will therma li ze it. We argue 
in Sec. [V] that for a certain experimentally relevant range of parameters, the thermalization 
is not significant over the relevant device length scales, justifying the approximation the 
injected electrons retain their 5-function distribution throughout the system. 

In the following sections we investigate the unstable modes which result from systems 
with an injected non-equilibrium beam of charged carriers. Before we present the results 
for graphene, we discuss the more familiar case of instability in three- and two-dimensional 
solid-state systems with parabolic energy band. We then compare the dispersion relation, 
region and strength of instability of the parabolic band and graphene cases. 

III. PARABOLIC BANDS 

In this section, we review the two-stream instability for parabolic bands where E(p ) = 
p 2 /2m, where m is the effective mass. The three-dimensional case has been well-studied 
by plasma physicists as it corresponds to a standard classical plasma where one charge 
distribution is streaming relative to another.- The two-dimensional case which corresponds 
to injection of a non-equilibrium beam of electrons into a two-dimensional electron gas such 
as a doped quantum well has also been considered by various authors.- - — ^ 

For parabolic bands, the polarizability of the equilibrium distribution subsystem 1 in the 
non-zero oj, small q limit is^S 



( 4 ) 


and, as shown in the appendix, the polarizability of the injected beam of charges is 



( 5 ) 
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where 9 is the angle between q and Vo, the velocity of the injected carriers. 

Using Eqs. flj]) and (jSJ) in Eq. ([[]) and setting e(q, u) equal to zero, gives the dispersion 
relation for three dimensional systems with parabolic band of 

47re 2 


1 = 


ni 

o 1 


n 2 


mix \_u] 2 (u — qvo cos 9) 2 


and for two dimensional systems of 

2vre 2 g 


1 = 


mix 


ni 


+ 


n 2 


( 6 ) 


( 7 ) 


; 2 (u — qvo cos#) 2 

where k is the dielectric constant of the background material. We note that these dispersion 
relations can also be obtained using Newton’s equations of motion and Poisson’s equation, 
which is the standard approach that is used in classical plasma physics.- The random-phase- 
approximation formalism that we have used, however, can be generalized to include quantum 
effects (which are important when q is not small), while the classical plasma physics method 
cannot. 

Comparing Eqs. (j6j) and (I7|) . there is an additional factor of q/2 on the right hand side 
of Eq. (jTjh which arises from the different forms of the Fourier transform of the Coulomb 
interaction in two and three dimensions. This factor of q/2 results in significant differences 
between the instabilities in three and two dimensions. 

In Eqs. (EJ) and (ITI) , the equilibrium value of the equilibrium charge density n\ is not 
necessarily equal to that of the injected charge carriers n 2 . Typically, in standard neutral 
plasmas, the densities are equal, since ri\ usually corresponds to the density of the positive 
ions and n 2 to the density of electrons. Since we are dealing with solid state systems, we 
assume that there is a uniform background charge that compensates for any net charge of 
the mobile charge carriers including injected charge. 

Multiplying Eqs. (j6j) and ((7j) by their respective denominators of the terms on the right 
hand side gives a fourth order polynomial equation for cv, so four roots are expected. The 
instability of the plasma wave is determined by whether the roots are all completely real 
or only two are completely real and two have non-zero imaginary parts. If only two roots 
are completely real, the remaining two roots appear as a complex conjugate pair. The root 
with a positive imaginary value of uj (he., Im(o;) > 0 ) indicates the exponentially growing 
wave, as can easily be seen by inserting a complex c o into the time-dependence of the wave, 
exp(— iojt), which results in a factor of exp[Im(o;)f]. Thus, Im(o;) gives the growth rate of 
the unstable wave.- 
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The solutions of the roots of Eqs. (J6j) and (J7J) were obtained numerically using Matlab. 
We assumed the ratio n\jn 2 = 10, where n\ = 10 18 cm -3 and ni = 10 12 cm -2 for the three- 
and two-dimensional cases, respectively. 

Fig. 1(a) and Fig. 2(a) show the imaginary part of the complex root of Eqs. © and 
(ED, respectively. It can be seen that for both three-and two-dimensional systems; the 
instabilities occur at low wavenumber and disappear above a certain cutoff q. In three- 
dimensional systems, the peak value of Im(u;) is the same for any given angle, but the peak 
wavenumber shifts to higher q as the angle 9 between q and po increases. The reason for 
this evident when the dispersion relation Eq. ((HD is examined. The only dependence on the 
dispersion on q is in the denominator of the second term on the right, where q is multiplied 
by cos 9. Therefore, the effect of changing the angle 9 is to “renormalize” the value of q 
for 9 = 0° to qcos9 for 9 7 ^ 0°. In the two dimensional case, in addition to a shift in the 
position of the peak of the growth rate Im(u;) towards increasing value of q with increasing 
angle 9 that is seen in the three-dimensional case, the magnitude of the peak also increases. 
This is due to the additional factor of q on the right-hand side of the dispersion relation. 
Figures (b) and (c) of both Figs. 1 and 2 show the real and imaginary parts of the dispersion 
relations for a representative small (10°) and moderately large (45°) angle of q with respect 
to po- There are four modes, two of which have frequencies that are always real, which 
correspond to the standard plasmon propagating in opposite directions, and two of which 
come in complex conjugate pairs below a certain q (which depends on the angle). 

These results for the growth rates are reproduced as surface and contour plots in Figs. 3 
and 4 for three dimensional and two dimensional parabolic band system respectively, where it 
is assumed that the momentum of the injected charge carriers is in the positive x-direct ion. 
We see easily from these plots that in the three-dimensional case the instability depends 
only on the component of q x (the component of q in the direction of charge carrier injection) 
and is independent of q y (which, by symmetry is equivalent to any component of q that is 
perpendicular to the x-direction). This is easily explained by inspection of the dispersion 
relation Eq. (|6j) , where the only dependence on the wavenumber q occurs as qcos9, which 
is equal to q x . In the two-dimensional case the growth rate is also dependent of q y , because 
of the presence of the additional factor of q — (g 2 + g 2 )^ in the dispersion relation Eq. (ED, 
as compared to Eq. (J6J) for the three-dimensional case. 

The plots presented in Figs. 3 and 4 also explain what appears to be a contradictory result 
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in Figs. 1(a) and 2(a). The range of wavevectors for unstable waves grows as 9 increases 
from 0° to 90° but appears to disappear abruptly and completely when 9 = 90°. Figs. 3 and 
4 show that there is no discontinuity when 9 —y 90°. What has happened in Fig. 1 and 2 for 
6 = 90° is that the peak of the instability has been pushed to q —> oo, therefore it does not 
appear on the plot. All that is visible is the q —> 0 dependence of Im(o;), which vanishes. 
We remind the reader that we have used the small q approximation and therefore the results 
presented here are not reliable at large q. 


IV. GRAPHENE 


We now investigate the instability when the energy dispersion is E(p) = pv 0 in a two- 
dimensional system as in the case of graphene, where the momentum is taken with respect 
to the K and K' points in the Brillouin zone. The velocity is taken to be Vq = 1.0 x 10 8 cm/s 
as in graphene. We take the density of the ri\ = 10 12 cnr 2 , which gives a Fermi energy of 
e F — hvoy/Ttrii = 120 meV, and the ratio of the density of the equilibrium electrons in the 
system to the injected electrons to be ri\ /n 2 = 10. We assume that the electrons are injected 
into the system at energy of E = lOmeV above the Fermi energy. 

For the equilibrium electrons at zero temperature, where the distribution function is 
constant for |p| < pp and zero for |p| > pp, the polarizability for the equilibrium electrons 
in the small q limit is^l 


rii(g,(u) = 


niVpg- 

PpOJ 2 


( 8 ) 


and, as shown in the appendix, the polarizability for a distribution of density n 2 that is 
peaked around po which makes an angle 9 with respect to q is 

n 2 Vpq 2 sin 2 9 


U 2 {q,uj) = 

Po((jU — qv o cos 9) z 

Thus setting the dielectric function equal to zero in Eq. (JT]) gives the dispersion relation 


(9) 


1 = 


27re 2 voq 


ni 1 

n T 


ri2 


sin 2 9 1 


( 10 ) 


k [uj 2 pp (u — qv 0 cos 9) 2 p 0 
where k is the dielectric constant in the 2-dimensional layer. Eq. no can be rewritten in 
the form 


1 


e 1 

+ 


£ 2 (z - A) 2 ’ 


( 11 ) 










where 


nipo 

n-ipF yin 2 9 


(12a) 


z 


= CO 


np 0 


27ie 2 Von2q sin 9 


CO 

co*(q ) 


(12b) 


and where 


A 


KPoVoq , a qv OCOS 9 

--— cot 9 =- 

27re 2 n 2 cu*(g) 


u\q) = 


'27re 2 v 0 ri2q 
np 0 


sin 9 , 


(12c) 


(12d) 


is the long wavelength plasmon dispersion for a density of carrier ri 2 that is injected into the 
system with momentum p 0 . This results in a quartic equation in z 


z 4 - 2\z 3 + (A 2 - e - 1 )z 2 + 2e\z - A 2 e = 0 . 


(13) 


Solving analytically, two of the four roots have non-zero imaginary parts when 2 - A < 
(1 + e 1 / 3 ) 2 / 3 . Since the coefficients of the quartic equations are real, the complex roots come 
in complex conjugate pairs. In Fig. 0 we plot the imaginary part of the complex root of 
Eq. (TT3j) (i.e. Im(z) as a function of A and e). In order to obtain the instability growth 
rates as a function of q = (q x , q y ) for certain material parameters, one first obtains A and e 
for that q and the material parameters. The Im[z(A, e)] is obtained from Eq. (ITT]) . and the 
growth rate can then be obtained from Eq. (I12b|) : i.e., Im(u;) = a/(g)Im(x) where co*(q) is 
given by Eq. (Il2dj) . The results for certain parameters are shown in Figs. 6 and 7. As in 
the case of the instability in the parabolic-band case in two dimensions, for a given q x , the 
growth rate increases with increasing q y , and the instability disappears for q perpendicular 
to the momentum of the injected carriers, Pq. However, there is an important qualitative 
difference between the parabolic band and graphene cases. For parabolic bands, instabilities 
occur (up to a certain wave-number) for q along the x-axis; that is, parallel to the direction 
of injection of the carriers. In the case of graphene, there is no instability along the x- 
axis. Mathematically, this is because of the presence of the sin 2 9 term in the polarizability 
n 2 (< 7 ,cu) for the injected carriers graphene, which is absent in the equivalent polarizability in 
the parabolic-band case. As in Figs.l (b), (c) and 2 (b), (c), the Fig. 6 (b) and (c) show the 
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real and imaginary parts of the dispersion relations for a small angle (10°) and a moderately 
large angle (45°) of q with respect to po respectively. The major difference between the 
graphene and the parabolic-band cases occurs for small angles. Note that in the graphene 
case, for the 10° graph, the scale of the wavevector and growth rates are roughly an order 
of magnitude smaller than for the 45° graph. In the next section, we discuss the reason for 
this significant qualitative difference between the parabolic-band and graphene cases. 


V. DISCUSSION AND CONCLUSION 

The absence of instabilities in modes with q in the same direction of the injected carriers 
Po in graphene is due to the linear energy dispersion. To understand why this is so, consider 
the mechanism for the formation of unstable modes. When there is a small spatial fluctuation 
in the charge density (caused by, for example, a thermal fluctuation) the resultant electric 
potential perturbation due to this fluctuation causes a change in the motion of the charge 
carriers. In systems that are in equilibrium, the net effect of the change in the motion of 
the charge carriers is to reduce the magnitude of the spatial fluctuation in charge. However, 
under certain non-equilibrium situations, as in the case of two-stream instabilities, a spatial 
fluctuation in the charge density causes a change in the motion of the charge carriers which 
tends to increase the fluctuation, resulting in a charge fluctuation that grows (initially) 
exponentially. Note that in order for this unstable feedback loop to occur, the spatial 
fluctuation in the charge density has to change the motion of the charge carriers. 

In the case of a one-dimensional system where the energy dispersion is linear, the carriers 
in an injected beam centered on momentum po, say in the positive direction, all have the 
same velocity v 0 . As the velocity of the particles is independent of the momentum, any 
change in the momentum of the particles due to the forces caused by the electric potential 
caused by a spatial fluctuation in the charge will not change the motion of the particles. 
Thus, the unstable feedback loop never occurs and there is no instability. 

This explains why there is no instability in two-dimensional graphene system when q 
is in the same direction of the injected carriers po, since this situation mirrors that of the 
one-dimensional system. But since the speed of particles is the same for any momentum, 
why does this argument not work when q is not in the same direction as p 0 ? As can be 
seen from the results shown in Figs. 6 and 7, there are robust regions of instability for q not 
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parallel to po- 

This is because while the speed of the particles is constant, their velocities (which also 
include the direction of motion) depend on the momentum. Therefore, when q is not parallel 
to p 0l the density perturbation of the wave produces a force which affects the direction of the 
particles in the beam centered at po, which affects the velocity of the particles. This allows 
the unstable feedback loop to occur, resulting in (initially) exponentially growing waves. 

We note that in this paper, we have made several simplifying assumptions. We have used 
a small wavevector approximation, which will break down when at wavevectors that are 
on the order of the smallest characteristic wavevectors of the distributions of the electrons. 
Therefore, the unbounded increase in the instability growth rates with increase in the mag¬ 
nitude of the wavevectors that are seen in the two-dimensional instabilities are expected to 
be an artifact of the small wavevector approximation. We have also not taken into account 
the effect of the lattice scattering on the distributions and the instability rates. The lattice 
scattering will tend to cause diffusive behavior in the small wavevector regime, which will 
tend to suppress the instability growth rates. 

Experimentally, one possible method for obtaining the non-equilibrium distribution of 
injected carriers in doped graphene that is required for the occurrence of the two-stream 
instability is to use a planar tunnel junction attached to an edge of a doped graphene 
flake.— The detection of the unstable waves could be achieved by placing, at the opposite 
side of the graphene flake from the injection end, several individually contacted drains which 
are at different angles with respect the direction of the injected electrons. There should be 
a variation in the magnitudes of currents detected by drains, with the largest enhancement 
occurring for gates that are placed at roughly 45° with respect to the direction of injection 
of the electrons, since the growth rates of the unstable waves in those directions are the 
largest. Since the frequencies of the unstable waves range up to the terahertz regime, these 
waves may be beyond the usual detection limit of standard electronic devices. In this case, 
it may be possible to detect the waves via a grating coupler in the drain region, which is 
designed to couple to the terahertz plasma waves and emit terahertz radiation, 24 which can 
then be detected.— In fact, using this technique, this instability might potentially be used 
as a source of terahertz electromagnetic waves. 

The possibility of observing two-stream instabilities in graphene are enhanced over stan¬ 
dard parabolic-band two-dimensional electron gases. In parabolic band systems, in order 
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to obtain the conditions for a two-stream instability, the non-eqnilibrium electrons must 
be injected at high energy to obtain the necessary speed. This makes them susceptible to 
the strong inelastic scattering processes that typically occur at high energies (such as optic- 
phonon and electron-electron scattering), which quickly degrade the current. In graphene, 
the speed of the electrons in is independent of the energy, and therefore there is no need to in¬ 
ject the carriers at high energies in order to obtain the conditions necessary for a two-stream 
instability. 

This fact supports our approximation that the non-equilibrium distribution function 
maintains its form (he., it is strongly peaked at the injection energy) as it passes through 
the system. The main contribution to the thermalization is the electron-electron scatter¬ 
ing of the injected carriers with the equilibrium carriers that are already in the sample. 
The single-particle lifetime of an electron in graphene due to electron-electron scattering 
with a fermi sea of electrons at zero temperature r ee is given by™ r” 1 = — 2Im[E ret ] = 


4 nep 


log 


8e_F 

s 


, where E ret is the on-shell retarded self-energy, £p is the Fermi energy 


and £ is the energy of the electron relative to ep. For £ = lOmeV and Sp — 120 meV, 
corresponding to the parameters used in Figs. (J6l) and (J7J), this gives r ee = 2.4 ps, which 
translates to a mean free path of 2.4 /irn. This is considerably larger than the wavelengths 
a good portion of the unstable modes. Furthermore, the r ee represents the single-particle 
lifetime, which takes into account all scattering out of a given momentum state, regardless 
of how small the momentum transfer is. Due to the long-ranged nature of the Coulomb 
scattering, the small-momentum scattering events dominate, but these events are very inef¬ 
ficient at degrading the non-equilibrium distribution function. A more appropriate measure 
of the rate of thermalization of the distribution is the transport lifetime, which weights 
scattering events by the momentum loss of the electron. We are not aware of any calcula¬ 
tions of the transport relaxation rate due to electron-electron scattering in graphene, but 
if we assume that charged-impurity scattering rates are an acceptable proxy for electron- 
electron scattering, then the transport lifetime can be over an order of magnitude larger 
than the single-particle lifetimes.— This furthers bolsters the validity of our approximation 
of a constant distribution function for the injected electrons across the device.— 

Finally, we contrast the two-stream instability described here with the more well-known 
Dyakonov-Shur (DS) instability,— which has been studied experimentally in field-effect 
transistors^ and theoretically in graphene.—The DS instability occurs in two-dimensional 
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electron gases, such as field effect transistors, which have a certain range of drift velocities 
and are subject to boundary conditions of constant total current at one end and constant the 
electric potential at the other. This leads under appropriate conditions to the amplitudes 
of density perturbations inside the two-dimensional electron gas being enhanced when they 
reflect off the boundaries, resulting in unstable waves. The two-stream instability on the 
other hand is a bulk effect, and the boundaries of the sample do not drive the effect. An¬ 
other major difference is that the DS instability is based on a hydrodynamic theory, which 
assumes that the electron-electron scattering dominates to the extent that the distribution 
function is locally a drifted thermalized distribution. This is the opposite limit from the 
two-stream instability, where it is crucial that the electron-electron scattering is at a low 
enough level to prevent thermalization of the non-equilibrium injected electrons. 

In conclusion, we have studied two-stream instabilities in two-dimensional systems with 
linear energy dispersions, such as in graphene, in which a beam of electrons is injected into a 
doped system. The range of wavevectors in which instabilities occur in the case of graphene is 
qualitatively different from that of conventional three- and two-dimensional parabolic-band 
systems. In particular, there is a complete absence of instabilities in waves with wavevectors 
that are parallel to the direction of injection of the non-equilibrium electrons. 
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Appendix: Derivations of Eqs. ([5|) and ([9]) 

The expressions for the polarizability in Eqs. (J5]) and ([9]) are derived from substituting 
the sharply peaked distribution n(p) = n 2 5( p — po) into Eq. (j2J). Let q = gx, and p 0 make 
and angle of 6 with respect to the x-axis, so that po )X = po cos 6 and |po,_l| = Po sin#, where 
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p_L are the components of p that are perpendicular to the x-axis. Then, 


n 2 (q,w) = n 2 


dp x 


dp± 


d 

q dp~ 1 2 3 4 S 6 ^ Px ~ Po ’ x ' _ p °a) 

qv x (p x , p±) -<jj 


= n 2 


= n 2 q 


= n 2 q 


dp. 


J all p_|_ 

Q~r\ Px Px, o) 

OPx _ 

qV x (Px,P0,y ) - ^ 

{h) S(v -V ' dV ^Px,PO,y) 

&px 0\Px Po,x 


{qv x ( p 0 ) 


dpi 
dvx(Px, Po.i 
dp x 


(qv x (p x ,Po,y) -ujf 


(A.l) 


- Px — PO,x 


where in the above we have integrated by parts to obtain the third equality. 

Qv x 1 _____ 

For parabolic bands, v x = p x /m , so that = —• Substituting this into Eq. (IA.1I) 

op x m 

gives Eq. (15j) . For a linear bands in 2 and 3 dimensions, v x = Vq — = , and therefore 

VpI + pI 

rp^ 

= Vq- —p. Substituting this into Eq. (lA.ljl . and using |p 0 ,j_| = posing, gives Eq. (J9]). 

dv x 

Note that for linear bands in one-dimension, v x is constant and therefore —= 0. (This 

OPx 

reflects the fact that in one-dimensional systems with a linear band, the motion of the 
particles does not depend on forces it experiences.) Therefore, two stream instabilities do 
not occur in one-dimensional systems with linear bands. 


dvx 

dp x 
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Figures 



FIG. 1: (a) Imaginary part of oj (the growth rate) as a function of wave number for different angles 
6 between wavevector q and injected carrier momentum po, for a three-dimensional parabolic band 
system. The mass is taken to be the bare electron mass, k = 1, n\ = 10 18 cm -3 , riijn-i = 10 and 
the energy of the injected electrons is E = lOrneV above the Fermi energy. Figs, (b) and (c) show 
the dispersion relation for angles 9 = 10° and 8 = 45° respectively. The solid lines are the real part 
of |w| and the dashed lines are the imaginary part of the modes that split at around 1.0 x 10 6 cm -1 
and 1.5 x 10 6 cnr 1 in (b) and (c), respectively. 
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q (10 5 cm' 1 ) 


FIG. 2: Same as in Fig. 1 for a two-dimensional parabolic-band system with n\ = 10 12 cm 2 , and 
other parameters unchanged. 



FIG. 3: Growth rates of the unstable plasmons (Im[o;]) as a function of q for a three-dimensional 
parabolic-band system, shown as (a) a surface plot and (b) a contour plot. The parameters used 
are the same as Fig. 1. The dark flat region in (a), and correspondingly the region to the right of 
q x = 1.0 x 10 6 cnh 1 in (b), is where Im(cu) = 0; i.e., stable waves. In (b), the lowest contour lines 
are at 2 x 10 12 s -1 and the difference between successive contour lines is 2 x 10 12 s -1 . 
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FIG. 4: Same as in Fig. 3 for a two-dimensional parabolic-band system with n\ = 10 12 cnr 2 , and 
other parameters unchanged. In (b), the lowest contour lines are at 1 x 10 12 s _1 and the difference 
between successive contour lines is 1 x 10 12 s -1 . 



FIG. 5: (a) Surface plot of Im( 2 ) as a function of mathematical parameters A and e, where z is the 
complex root of Eq.(14); (b) corresponds the counter plot of (a), where the lowest contour line is 
at 0.1 and the difference between successive contour lines is 0.1. 
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q(10 6 cm" 1 ) 


FIG. 6: (a) Imaginary part of u> (the growth rate) as a function of wave number for different 
angles 9 between wavevector q and injected carrier momentum po, for a linear-band system such 
as graphene. Here, vq = 1.0 x 10 8 cm/s, k = 3 (as for graphene on a BN substrate), n\ = 
1.0 x 10 12 cm” 2 and ni/ri 2 = 10 and the energy of the injected electrons is E = lOrneV above the 
Fermi energy. For both 0 = 0° and 90°, the imaginary part of u is zero. Figs, (b) and (c) show the 
dispersion relations for angle 9 = 10° and 9 = 45° respectively; as in Figs. 1 and 2, the solid lines 
are the real part and the dashed lines are the imaginary part of the mode that splits. 
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FIG. 7: (a) Surface plot and (b) contour plot for Im(w) as a function of q for graphene, with the 
same parameters as in Fig. 6. In (b), the lowest contour line is at 5 x 10 12 s^ 1 and the difference 
between successive contour lines is 5 x 10 12 s _1 . 
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